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Charge conjugation (C) and Parity (P) are exact symmetries at ^ = tt and = p,/ (iT) — n, where 9 is the 
parameter of the so-called ^-vacuum, jj, is the imaginary quark-number chemical potential and T is the temper- 
ature. Spontaneous breakings of these discrete symmetries are investigated by the the Polyakov-loop extended 
Nambu-Jona-Lasinio (PNJL) model. At zero T, P symmetry is spontaneously broken while C symmetry is 
conserved. As T increases, P symmetry is restored just after C symmetry is spontaneously broken, so that either 
P or C symmetry or both the symmetries are spontaneously broken for any T. The chiral-symmetry restora- 
tion and the deconfinement transition at = 6> = are remnants of the P restoration and the C breaking at 
6 = — TT, respectively. 
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I. INTRODUCTION 

Violations of parity (P), charge conjugation (C) and charge- 
parity symmetries (CP) in strong interaction are one of impor- 
tant subjects in particle and nuclear physics. For example, the 
strong CP problem is a long-standing puzzle; see for example 
Ref. ild for a review of this problem. Lorentz and gauge in- 
variance allow the Quantum Chromodynamics (QCD) action 
to have a term 



(1) 



of the topological charge, where F^^^ is the field strength of 
gluon. The parameter can take any arbitrary value between 
— TT and TT, where 6 = — tt is identical with 9 = n. Never- 
theless, experiment indicates |0| < 3 x 10~^° Since 9 is 
P-odd (CP-odd), P (CP) is then preserved for 6* and ±tt, 
but explicitly broken for other 6. Why is 9 so small ? This is 
the so-called strong CP problem. 

For zero temperature (T) and zero quark-chemical po- 
tential (/i), P is conserved at = 0, as Vafa and Witten 
showed Jsj]. Meanwhile, P is spontaneously broken 9 = n, 
as Dashen and Witten [5] pointed out. This is the so- 
called Dashen phenomena. Since the spontaneous P violation 
is a nonperturbative phenomenon, the phenomenon was so far 
studied mainly with the effective model such as the chiral per- 
turbation theory ll6l-[Tlll. 

For T higher than the QCD scale Aqc d, there is a possibil- 
ity that a finite 9, depending on spacetime coordinates {t, x), 
is effectively induced llT2ll . since sphalerons are so activated 
as to jump over the potential barrier between the different de- 
generate ground states |13]. If so, P symmetry can be vio- 
lated locally in high-energy heavy ion collisions. This effec- 
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tive 9{t, x) deviates the total number of particles plus antipar- 
ticles with right-handed helicity from that with left-handed 
helicity. The magnetic field, formed in the early stage of 
heavy-ion collision, will lift the degeneracy in spin depending 
on the charge of particle. As a consequence of this fact, an 
electromagnetic current is generated along the magnetic field, 
since particles with right-handed helicity moves opposite to 
antiparticles with right-handed helicity. This is the so-called 
chiral magnetic effect (CME) lfl4l - [l6ll . CME may explain the 
charge separations observed in the recent STAR results IitIi . 
Thus, theoretical study on the thermal system with nonzero 9 
is interesting. 

For finite fi, the QCD action has a term 



r ^ - 
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of the baryon-number charge, where q is the quark field. 
When fi is pure imaginary, i.e. /i = iOT, has a mathemat- 
ical structure similar to Ce, if the baryon number is conserved. 
The dimensionless chemical potential O can vary from — tt to 
TT, where O = — tt is identical with O = n. Since 6* is a 
C-odd quantity, C is an exact symmetry at = and ±7r, 
but not at other 0. Thus, C violation induced by finite is 
analogous in principle to P violation induced by finite 9. 

At imaginary /i, QCD has a periodicity of 27r/3 in 0. This 
periodicity was found by Roberge and Weiss with perturba- 
tive QCD for high T and strong-coupling QCD flSj for low 
T. This Roberge- Weiss (RW) periodicity was confirmed by 
lattice QCD (LQCD) lfT4j28ll . At higher temperature, there 
exist three Z3 vacua. As increases from — tt to tt, the three 
vacua emerge one by one. As a consequence of this mecha- 
nism, three first-order phase-transitions appear at 6* = ±7r/3 
and TT. At 6* = ±7r/3 and tt, thus, a mechanism similar to the 
Dashen phenomena at 6* = tt takes place. The transitions at 
= ±7r/3 and tt are called the RW transition. The C break- 
ing at 6* = ±7r/3 and tt occurs when T is high [29J, while 
the P breaking at 9 = n takes place when T is zero lUt] and 
then small. The grand canonical partition function at 6* = 
is a sum of the canonical partition function over the quark 
number that is obtained by the Fourier transform of the grand 
canonical partition function with finite U8ii30il . This means 
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that the singular behavior of the RW transition in the vicin- 
ity of 6* = 7r/3 reflects on the behavior of QCD at 6* = 0. 
Actually, it is confirmed in Ref. 1E9I1 that the deconfinement 
crossover at = is a remnant of the first-order RW transi- 
tion at 6* = 7r/3. Thus, the thermodynamics at = 7r/3 is 
closely related to that at 6* = 0. Furthermore, we can expect 
from the analogy between Cg and £^ that the thermodynam- 
ics at 6* = TT is also closely related to that at 6* = 0. 

LQCD has the sign problem at finite 9, but not at finite 0. 
Therefore, we can test an effective model at 6* = and 6* > 
and apply the model to the case of 6* > and 63 > 0. As a can- 
didate of such effective models, we can consider the Nambu- 
Jona-Lasinio (NJL) model Il3ll - l4lll and the Polyakov-loo p ex - 
tended Nambu-Jona-Lasinio (PNJL) model ifm l29l l42l468ll . 
The NJL model can describe the chiral-symmetry breaking 
and the Dashen mechanism ll39l - l4Til . but not the confinement 
mechanism and the RW transition. The PNJL model can treat 
the deconfinement and the RW transition |29, 54, 55] as well 
as the chiral symmetry breaking and CME [16] . However, 
the PNJL model has a weak correlation (entanglement) be- 
tween the chiral and the deconfinement transition compared 
with LQCD |55, 65]. In order to solve this problem, we re- 
cently proposed a new version of the PNJL model, i.e. the 
entanglement-PNJL (EPNJL) model {6^, that has a four- 
quark vertex depending on the Polyakov loop. The EPNJL 
model can reproduce not only the strong correlation between 
the two transitions without L25> .66,] and with the strong mag- 
net field ll68l I69I1 but also the quark-mass dependence of the 
order of the RW endpoint HHM]. 

The chiral and the Z3 symmetries are not exact symmetry 
for QCD with physical quark mass. Hence, the chiral and de- 
confinement transitions can not be defined exactly. Indeed, 
these transitions are approximately defined by the chiral con- 
densate and the Polykov loop. In addition, the two transitions 
are crossover at 6* = = 0. These situations often make 
it complicated the relation between the two transitions. At 
6 = = TT, in contrast, P and C are exact symmetry, so that 
their spontaneous breakings are clearly defined by P-odd and 
C-odd quantities respectively. Thus, the region of = 6* = tt 
is suitable to investigate the interplay between two kinds of 
transitions. Furthermore, we can expect that the thermody- 
namics at 6* = 6* = TT is closely related to that at 6* = 69 = 0. 

In this paper, we analyze P and C violations at finite 9 and 
and the interplay between them, using the PNJL and EP- 
NJL models. Particularly at 6* = 6? = tt, P and C are exact 
symmetries, so the analysis is mainly focused on the region. 
We also investigate the relation between the P and C breakings 
?X9 = = t: and the chiral and deconfinement transitions at 
9 = = Q. The P breaking at finite 9 was already stud- 
ied for the case of T = by the NJL model |39-41], so the 
present analysis is concentrated on the case of finite T, since 
the PNJL and the EPNJL model ai-e reduced to the NJL model 
in the Hmit of T = 0. 

This paper is organized as follows. In section II, the PNJL 
model is explained briefly. In section III, C and P violations 
are numerically investigated, particularly at 6? = 7r/3 and 9 = 
it; note that 6? = 7r/3 is identical with = n because of the 
RW periodicity. We also explain the EPNJL model briefly in 



this section. Section IV is devoted to summary. 



II. PNJL MODEL 

Pioneering work on the parity violation and its restoration 
in the framework of the NJL model was done by Fujihara, In- 
agaki and Kimura J39J. Boer and Boomsma studied on this 
issue extensively 1401 14111 . Here, we extend their formalism 
based on the NJL model to that on the PNJL model. The 
two-flavor {Nf = 2) PNJL Lagrangian with the 6'-dependent 
anomaly term is given as 

C = q{i^,D- -ni)q-U{<P{Al<P{A\\T) 



a=0 

+ 8G2 [e'' det [qnqL) + e'^' det (qi^qn)] , (3) 

where q = {u, d) denotes the two-flavor quark field, m does 
the current quark-mass matrix diag(mu, m^), tq is the 2x2 
unit matrix, Ti{a = 1,2,3) is the Pauli matrices and D'^ = 
d" + iA" ~inS^. The field A" is defined as A" = S^gA^^^ 
with the gauge field the Gell-Mann matrix Aa and the 
gauge coupling g. In the NJL sector, Gi denotes the cou- 
pling constant of the scalar and pseudoscalar-type four-quark 
interaction, and G2 is the coupling constant of the Kobayashi- 
Maskawa-'t Hooft determinant interaction ItoIItiIi the matrix 
indices of which run in the flavor space. The Polyakov poten- 
tial U, defined later in ( |20| ). is a function of the Polyakov loop 
<P and its Hermitian conjugate <P*, 



1 



TrX, <5* 



1 



TrL^ 



with 



r /-/^ 

L(x) =Vexp i / drA4(x, r) 
L Jo 



(4) 



(5) 



where V is the path ordering and A4 ~ iAg. In the chiral 
limit (TOu = = 0), the Lagrangian density has the exact 
SU{Nf )i^ X SU{Nf )n x U{1)^ x SU{3)c symmetry. The 
C/(1)a symmetry is explicitly broken if G2 7^ 0. The tem- 
poral component of the gauge field is diagonal in the flavor 
space, because the color and the flavor space are completely 
separated out in the present case. In the Polyakov gauge, L 
can be written in a diagonal form in the color space 14411 : 

^ ^ gi/3(03A3+08A8) ^ diag(e*'''^",e*''*%e*'^'^'^), (6) 

where 0a = 03 + ^s/VS, (t>b = -03 + ^s/VS and 0c 
— (0Q + 0h) = — 208/\/3. The Polyakov loop <P is an exact 
order parameter of the spontaneous Z3 symmetry breaking in 
the pure gauge theory. Although the Z3 symmetry is not an 
exact one in the system with dynamical quarks, it may be a 
good indicator of the deconfinement phase transition. There- 
fore, we use <P to define the deconfinement phase transition. 

For simplicity, we assume below that rriu = nid = mo. 
Furthermore, to remove the 9 dependence of the determinant 
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interaction, we transform the quark field q into the new field 

q' as 



(7) 



Under this U(1)a transformation, the quark and antiquark 
condensates are also transformed as 



= cos [ - I cr + sm I - ) 77 



9*759 = ^^^^12/'^ + cos I - ) 



qnq = cos [ - I a, + sm 12)^" 



qinjsq = - sm [ - l a, + cos l - ) tt, 



(8) 



where a' is defined by the same form as a but q is replaced 
by q'; this is the case also for other condensates rj', a'^ and tt'. 
The Lagrangian density is rewritten with the new field q' as 

£ = ct{i^,D^ -ma+-mM,^r^)q' -U{$[AmA]\T) 



Gi [{q'raq'r + {q'nr.raq' f] 



a=0 



(9) 



+ 8G2[det(g\90+dct((z\gk)] 

where 7710+ = 777,ocos(|) and 7719- = 7779 sin (|). Making 
the mean field approximation and performing the path inte- 
gral over the quark field, one can obtain the thermodynamic 
potential H (per volume) for finite T and fi: 

d3p r 



J? 



(2^)3 



3{E+{p) + E^ip)} 



In [1 + 3<Pe 



-pE-_ 



3^>*e 



-2/3 _b: 



iln [l + M*e-^''t +Me-^^''t 

iln [l + 3^*e-^^- 
U+U, 



-3/3B+] 



3^e 



where = E+ (p) ± fi and = S_ (p) ± ^ with 



E± - Vp2+C±2\/Z), 



(10) 

(11) 
(12) 



D = A^M^ + P^N^ + 2APMN cos ip + A^P^ sin^ LP 



{MA + NPf + (A X P)2 > 



(13) 



M ^ mo+ - 2G+r7' = 7770+ - 2(Gi + G2)a', (14) 

iV = 7770- - 2G_77' = 7770- - 2(Gi - G2)??', (15) 

A (-2G_a;,-2G_4,-2G_a3), (16) 

P - (-2G+7r;,-2G+7r^,-2G+7r^), (17) 

A = VA-A, P = VP • P, A-P = AFcos^, (18) 

U = G+ia''+n'^') + G^ia'^'+rj''). (19) 



In the right-hand side of (fTOl i. only the first term diverges. The 
term is then regularized by the three-dimensional momentum 
cutoff A HEa]. Following Ref. Hlill], we introduce c as 
Gi = (1 - c)G and G2 = cG, where < c < 1 and G > 0. 
Hence, the NJL sector has four parameter of ttio, A, G and 
c. We put 7770 = 5.5MeV. The parameters A and G are so 
chosen as to reproduce the pion decay constant /jr = 93MeV 
and the pion mass 7777r = 139MeV at vacuum. The remaining 
parameter c is a free parameter. Although the exact value of c 
is unknown, it is known from the analysis of the 77-77' splitting 
in the three flavor model that c ^ 0.2 is favorable |72]. The 
value c — 0.2 has been also used in Refs. I40il4l(l . Therefore, 
we adopt c = 0.2 in this paper (Note that the PNJL and the 
EPNJL model are reduced to the NJL model in the limit of 
T = as is mentioned above. ) For comparison, we will also 
show the result of the NJL model that has the same parameter 
set as in the PNJL model. 

The Polyakov potential U of Ref. jilis fitted to LQCD 
data in the pure gauge theory at finite T ll73ll74ll : 



U = T* 



a{T) 



b{T) ln(l - e-Zx?* + 4(<Z>=^ + <Z>*^) - 3(^'<Z' 



a(T) = ao + ai(^)+a2(^)^ h{T)^b,[^ 



(20) 
(21) 



The parameters included in U are summarized in Table I. The 
Polyakov potential yields a first-order deconfinement phase 
transition at T = To in the pure gauge theory. The original 
value of To is 270 MeV evaluated by the pure gauge lattice 
QCD calculation. However, the PNJL model with this value of 
To yields somewhat larger value of the transition temperature 
at zero chemical potential than the full LQCD simulation iItsI - 
FttIi predicts. Therefore, we rescale Tq to 212 MeV in the 
numerical calculations |d2|1. 



ao 




02 


bi 


3.51 


-2.47 


15.2 


-1.75 



TABLE I: Summary of the parameter set in the Polyakov-potential 
sector used in Ref. 14811 . All parameters are dimensionless. 



The variables X = <P, (p* and a satisfy the stationary con- 
ditions. 



df2/dX = 0. 



(22) 



The solutions of the stationary conditions do not give the 
global minimum f2 necessarily. There is a possibility that they 
yield a local minimum or even a maximum. We then have 
checked that the solutions yield the global minimum when the 
solutions X{T, 9, 0) are inserted into ( fTOl i. 

Now we consider the imaginary chemical potential ^ — 
iOT. The thermodynamic potential Q of ( fTOl i has a period- 
icity of 27r in both 6 and 0. Hereafter, we mainly consider 
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one circle, < 61 < 27r and < 6* < 2tt. In addition, f] has 
the RW periodicity: 

n{T, e, 0) = n{T, e,0 + ^) = q{t, 0,0 + ^) (23) 

for — 27r/3 < 6* < 0. This is understood as follows. The ther- 
modynamical potential Q is not invariant under the Z3 trans- 
formation, 



(24) 



while U of (|20] | is invariant. Instead of the Z3 symmetry, 
however, Q is invariant under the extended Z3 transforma- 
tion l.55i] . 



±10 



(25) 



for 6* = and tt. The potential J? is also invariant under C 
transformation. 



7r2 



-7r2, 02 



-02, 



(29) 



for = and tt. Because of the RW periodicity (|23i l, the C- 
invariance is true also for = mr/i, where n is an arbitrary 
integer Hereafter, we mainly consider a period < < tt/3 
for simplicity. 

For — 7r/3, C-symmetry is spontaneously broken at T 
higher than some critical temperature Tc ll29ll . C-odd quan- 
tities such as the phase of if", the imaginary part of ^P or the 
baryon number density is the order parameter of this phase 
transition. For 6 ^ tt, as shown later, P-symmetry is spon- 
taneously broken when c is greater than some critical value 
Ccri and T is smaller than some critical temperature Tp JsH- 
I4TI1 . P-odd quantities such as 77 are the order parameter of this 
phase transition. 



This invariance means that f2 has the Roberge-Weiss period- 
icity 05l]. This can be seen more explicitly by introducing the 
modified Polyakov loop 



^ = (Pe 



i0 



(26) 



invariant under the extended Z3 transformation dZST l. The ther- 
modynamic potential Q is then rewritten with the modified 
Polykov loop as 



(27r)3 



3{£;+(p) + i?_(p)} 



3^*g-2/3£;+(p)+i3e _^ g-3/3£;+(p)+3' 



101 



3^*g-2/3£;_(p)-3ie _^ g-3/3£;_(p)-3iej 



+ 3i^/e-2/3S+(p)+3je ^ g-3/3B+(p)+3i0j 

+ ihi [l + 3tf'*e-'^^-(P) 

P 

+ 3iZ/e-2/3S-(p)+3»e g-3/3£;_(p)+3ie] 



a{T) 



W*W + b{T) ln(l - 6if If * 



3„-3j0 



*\2^ 



In (|27] |. depends on only through e*''®. Thus, has the 
RW periodicity (l23T l. 

The thermodynamic potential ]7 is invariant under P trans- 
formation. 



(28) 



III. P AND C BREAKING AT FINITE 9 AND 6> 

In this section, P and C violations are investigated first by 
the PNJL model. After the PNJL analysis, we use the EPNJL 
model to explain the strong correlation between the chiral and 
deconfinement transitions at zero and finite 6 and 0. 

The present PNJL model has eight condensates of quark- 
antiquark pair. However, a and 7? vanish ll40ll4Tll . since niu = 
rrid and the isospin chemical potential is not considered here. 
We can then concentrate ourselves on a, r/ and <P. 



A. Thecaseof = 6> = 

In this subsection, we consider the case of 6 — — 0. 
Figure [TJa) shows T dependence of the chiral condensate a 
and the Polyakov loop P; here, a is normalized by the value 
ctq at T = 0. Rapid but smooth changes in a (solid line) and 
P (dashed line) indicate that the chiral-symmetry restoration 
and the deconfinement transition are crossover Figure [Ttb) 
presents the chiral and Polyakov-loop susceptibilities, Xa-a 
and x<fs>*, as a function of T. The pseudocritical tempera- 
tures, and Td, of the chiral-symmetry restoration and the 
deconfinement transition are defined by peaks of Xa-a- (solid 
line) and x*** (dashed line); here, the susceptibilities are nor- 
malized by T so as to be dimensionless p4]. In the present 
PNJL model, = 216 MeV and Id. = 173 MeV, while 
n « Td 



173 MeV in LQCD U3\M- The value of is 



even larger than the NJL result T^JL = igg MgV. Thus, the 
present PNJL result is consistent with LQCD data for Td, but 
(27) not for T^. This will be discussed at the end of this section. 



B. The case of 6» > and = 

In this subsection, is fixed at zero. Figure |2 a) shows T 
dependence of cr, (p and the absolute value of 77 at 6* = tt; here, 
a and |?7| are normalized by |ryo| the value of |r/| at T = 0. 
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Fig. 1: (color online). T dependence of (a) the chiral condensate 
and the Polyakov loop and also (b) their susceptibilities at 6* = 
and = 0. The solid (dashed) line represents the chiral condensate 
(Polyakov loop) in (a) and the chiral (Polyakov-loop) susceptibility 
in (b). In panel (b), x*** is multiplied by 10'^. 



The solid and dashed lines represent |?7/?7o| and respec- 
tively, while the dot-dashed line stands for a/ 1 ryo | ■ The order 
parameter 77 of P violation is finite at T < Tp = 202 MeV, 
while zero at T > Tp. This means that P-symmetry is sponta- 
neously broken below Tp and restored above Tp. The restora- 
tion of P breaking at T = Tp is a second-order transition, as 
shown by the solid line. T dependence of <!> little changes be- 
tween Fig.ina) and Fig. IS a), indicating that Td = 173 MeV 
also for 6 = TT. Thus, Tp is much higher than Td. Com- 
paring the solid and the dot-dashed lines shows that \a\ is 
much smaller than \ri\ at small T. Thus, the thermodynam- 
ics at = TT is mainly controlled by I77I. 

The previous works on the NJL model showed for 9 — tt 
that P symmetry can be violated for low T, but not for high 
T Il39l - l4lll . This statement is supported by the present PNJL 
model. In Fig.|2ja), the dotted line is a result of the NJL model 
for |?7/?7o|- The order of P restoration is second order [39-41 J, 
as shown by the dotted line. Thus, the PNJL and the NJL 
model show that the P restoration is a second-order transition, 
while the linear sigma model points out that it is a first-order 
transition IiIItI. In the present NJL model, Tp « 172 MeV 
and then close to Td, but this is just accidental. 

Figure|2|b) shows T dependence of S = y^a^ + rf at sev- 
eral values of 9. For 9 = 0, S agrees with \a\ the approximate 
order parameter of the chiral symmetry. For 6* = tt, mean- 
while, S is close to the order parameter jryj of P violation. 
Although S itself is not an order parameter of P restoration, 
it has a cusp at T = Tp as a reflection of the second-order 



P restoration at T = Tp. This nonanalytic behavior in E is 
smeared out as 9 decreases from tt, but a smooth but rapid T 
dependence remains in i7 at 6* = 0. Thus, the crossover chiral 
symmetry restoration at 9 = can be regarded as a remnant 
of the second-order P restoration ax9 = n. 
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Fig. 2: (color online), (a) T dependence of a, $ and the absolute 
value of rj at 6 — TV and also (b) a/ct^ + 77^ in T-6 plane. In panel (a), 
the solid (dashed) line represents 1 77/770 1 i'^), while the dot-dashed 
line stands for a/\rio\. The dotted line shows 177/770 1 calculated with 
the NJL model. In panel (b), the solid lines show the PNJL results at 
several values of 9. In both panels, = 0. 

Figure[3la) show 9 dependence of the thermodynamical po- 
tential n. We see that f2 is minimum at 9 = and is max- 
imum at 9 = n. When T < Tp (solid line), due to the P 
violation at 9 = tt, f2 has a cusp there. 

Figure |3jb) show 9 dependence of cr (solid and bold solid 
lines) and 77 (dashed and bold dashed lines). We see that, at 
low temperature, due to the P violation, a has a cusp at9 = n, 
while 77 is discontinuous there. In general, below the critical 
temperature Tp, 6'-even quantities such as i7 or cr have a cusp 
at = TT, while 6'-odd quantities such as rj are discontinuous 
there. 

FigureOc) shows the phase diagram in T-9 plane at 6? = 0. 
The vertical solid line represents the first-order phase transi- 
tion induced by spontaneous P breaking, since on the line J? 
is not smooth in the 9 direction. The dashed line stands for the 
chiral crossover connected with the solid line at the endpoint 
of P violation. The order of the endpoint is the second order, 
as mentioned above. The dotted line is the deconfinement- 
crossover line that is almost independent of 9. 
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Fig. 3: (color online), (a) 6 dependence of the thermodynamic po- 
tential. The dashed (solid) line represents the result at T = 220MeV 
(leOMeV). no = n{e ^ 0,0 = 0,T). (b) e dependence of a 
and rj. The solid (bold solid) line represents a at T = 220MeV 
(160MeV), while the dashed (bold dashed) line represents the 77 
at T = 220MeV (160MeV). (c) Phase diagram in T-O plane at 
= 0. The vertical solid line represents the first-order phase transi- 
tion caused by P violation. The dashed and dotted lines stand for the 
chiral and the deconfinement crossover, respectively. 



C. The case of > and 6* = 

In this subsection, we review the case of finite to see the 
analogy between P violation induced by finite 6 and C viola- 
tion induced by finite 0. For this purpose, 9 is fixed at here. 
Figure |4f a) shows T dependence of a, the absolute value of 
the phase tp of 'F, the imaginary part of the quark number 
density Uq a.t0 ~ 7r/3. As shown by the dashed (dotted) line, 
tjj (Im(ng)) is zero below Tc — 189 MeV but finite above 
Tc- This indicates that C symmetry is spontaneously broken 
above Tq |29]. The order of the C violation (RW transition) 
is of first order, since i7 is not smooth in the direction |55]. 
The order of the endpoint at T = Tc of the C violation line 



depends on the Polyakov potential U used. It is the second 
order jS for U of Ref. |44], but the first order L61J for U of 
Ref. l48ll . The agreement of the PNJL result with the LQCD 
data at finite is better in the later than in the former ll6l[l . 
Therefore, we take the latter case in this paper Hence, the 
endpoint of the C violation (RW transition) is a triple point 
in the present case; this is explicitly shown in Ref. f6V\. Be- 
cause of this first-order phase transition, cr (solid line) has a 
gap at T = Tc. This behavior of a is easily understood by 
the discontinuity theorem on the first-order phase transition 
by Barducci, Casalbuoni, Pettini and Gatto |79]. However, 
the absolute value of cr is still large at T = Tc, as shown by 
the solid line. Hence the chiral transition should be regarded 
as not of first order but crossover. Comparing Fig. |3a) with 
Fig. we can see that the chiral crossover restoration is 
slower at 6? = 7r/3 than at 6* = 0. 

Figure Ub) shows the absolute value of the Polyakov loop 
<P in T-0 plane; note that = by definition. The abso- 
lute value of <P is not an order parameter of C violation, but 
\<l>\ (solid line) has a jump near T = Tc and 6* = 7r/3 as a 
reflection of the fact that the endpoint of C violation is of first 
order at T = Tc and — tt/3. This nonanalytic behavior 
in 1^1 is smeared out as decreases from 7r/3. However, a 
smooth but rapid T dependence remains in \<P\. Thus, the de- 
confinement crossover at 6* = is a remnant of C violation at 
= ^/3 129]. 

Figure 13 a) show the dependence of the thermodynam- 
ical potential f2. The RW periodicity ( |23] | is clearly seen in 
this figure. We also see that i7 is minimum at 6* = (27r/3, 
47r/3 ) and is maximum value at 69 = 7r/3 (tt, 57r/3). When 
T > Tc (dashed line), due to the C violation at 6* = 7r/3 (tt, 
57r/3), SI has a cusp there. 

Figure |5jb) show the dependence of cr (solid and bold 
solid lines) and the imaginary part of the quark number den- 
sity Uq (dashed and bold dashed lines). We see that, at high 
temperature, due to the C violation, a has a cusp at = 7r/3 
(tt, 57r/3), while Im(nq) is discontinuous there. In general, 
above the critical temperature Tc, 69-even quantities such as 
SI or a have a cusp at 6* = 7r/3 (tt, 57r/3), while 6?-odd quan- 
tities such as Im(7ig) are discontinuous there. 

The phase diagram in the T-0 plane is shown in Fig.|5jc). 
The vertical solid line represents the first-order RW transition 
induced by spontaneous C-breaking, since on the line Q is 
not smooth in the direction. The dashed line stands for 
the chiral crossover The dotted line is the deconfinement- 
crossover line connected with the solid line at the endpoint 
of C violation. The order of the endpoint is of first order, as 
mentioned above. 



D. The case of S > and 6» > 

In this subsection, we consider the case of finite 6 and 
to examine the correlation between C violation and P restora- 
tion. Figure |6t a) shows T dependence of \rj\ at 6* = and 
7r/3 in the case of 6 ~ n. The solid and the dashed line show 
the PNJL results at = and 7r/3, respectively. The crit- 
ical temperature Tp of P-restoration is 202 MeV for = 
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Fig. 4: (color online), (a) T dependence of a (solid line), the absolute 
value of the phase if) of (dashed line) and the absolute value of 
the imaginary part of the quark nunber density Hq divided by NjT^ 
(dotted line) at 6* = 7r/3 and also (b) the absolute value of $ in T-0 
plane. In both cases, 6 = Q. 



and 249 MeV for O — tt/S. Thus, the imaginary chemical 
potential enhances the P violation as well as the chiral sym- 
metry breaking. Figure |6jb) shows T dependence of jr;! and 
at 6* = 7r/3 and 9 = n; here, \t]\ is normalized by the 
value 1 770 1 at r = 0. The solid and the dashed line stand 
for 1 77 1 and respectively. A first-order C violation occurs 
at Tc — 190 MeV, while a second-order P-restoration takes 
place at Tp = 249 MeV. Thus, Tp is much larger than Tq. 
This indicates that either P or C symmetry or both are broken 
for any T. The inequality Tc < Tp at 6* = tt and 6> = 7r/3 
shown in Fig.|6lb) and the inequality Td ^ T^ at 6* = and 
= shown in Fig. [Tfb) come from the fact that the cor- 
relation between the chiral and deconfinement transitions are 
weak in the present PNJL model. This problem will be solved 
in the next subsection by using the EPNJL model. 

The order of C violation depends on the Polyakov poten- 
tial U taken. As mentioned above, the order is of first order 
for U of ( |20b proposed by Rosner et al. ll48ll . Meanwhile, 
as shown in Fig. |71 it is of second order for U proposed by 
Fukushima L44.1 , 



U = -bT 



54e-"/^<l>*(Z> 



ln(l - 6<P<P* + 4:{<P^ + -P*^) - 3{<P<P*f) , (30) 
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Fig. 5: (color online), (a) O dependence of thermodynamic po- 
tential at ^ = 0. The solid (dashed) line represents the result for 
T =175(195)MeV. = J7(6l = 0, = 0, T). (b) dependence 
of (7 and Im(nq) at S = 0. The solid (bold solid) line represents ct at 
T = 175(195)MeV, while the dashed (bold dashed) line represents 
Im(7i,) at T = 175(195)MeV. (c)The phase diagram in T-O plane 
at 9 = 0. The vertical solid line represents the first-order phase tran- 
sition of C violation (the RW transition). The dashed line stands for 
the chiral crossover, while the dotted line shows the deconfinement 
crossover. 



where a = 664MeV and h = 0.015/1^ f2§\. However, the 
first order is more plausible, since lA proposed by Rosner et 
al. is more consistent with LQCD data at imaginary chemical 
potential than that by Fukushima 1I61I1 . Furthermore, 8-flavor 
LQCD data at imaginary isospion chemical potential |24] fa- 
vors lA proposed by Rosner et al, since the EPNJL model with 
lA proposed by Rossener et al. yields results consistent with 
the LQCD data 16^. 
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Fig. 6: (color online). T dependence of (a) at = and = 
7r/3 and also (b) |7?/?;o | and \'4>\ at = 7r/3. Here, 9 is fixed at vr. In 
panel (a), the solid (dashed) line represents the PNJL result at = 
(tt/S). In panel (b), the solid (dashed) line stands for rj (tp). 
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Fig. 7: (color online). T dependence of |?7/r;o| (solid line) and {tpl 
(dashed line) at ^ = tt and — 7r/3 in the PNJL model with the 
Polyakov potential proposed by Fukushima t441 . 



Z3 symmetry, 

G{<P) = G[l - ai^<P* - a2{<P^ 



*3\ 



(31) 



This model is called the entanglement PNJL (EPNJL). The 
EPNJL model with the parameter set, ai = a2 = 0.2 and 
To = 190 MeV, can reproduce LQCD data at imaginary 
chemical potential and real isospin chemical potential as well 
as the results at zero chemical potential [66J. The EPNJL 
model with this parameter set is applied for the present case 
with zero and finite 6 and 0. 

Figure [8] shows results of the EPNJL model. Panel (a) 
presents T dependence of a (solid line) and ^ (dashed Une) 
at 6* = 6* = 0, and panel (b) shows T dependence of Xaa 
(solid line) and x^<i>' (dashed line) at 6* = 63 = 0. The entan- 
glement vertex dSTT l makes the chiral and the deconfinement 
crossover almost coincide. 
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Fig. 8: (color online). Results of the EPNJL model for T dependence 
of (a) a I GO and ^ at 6* = = and (b) Xo-o- and x*** at ^ = = 
0. The solid (dashed) line represents ct (<?>) in panel (a) and Xo-ct 
(X***) in panel (b). In panel (b), x*** is multiplied by 10'^. 



E. The EPNJL model 

Since the four-quark coupling constant G contains effects 
of gluons, G may depend on In fact, recent calcula- 
tions 180148211 of the exact renormaUzation group equation 
(ERGE) [SJ] suggest that the chiral and the deconfinement 
transitions coincide by the higher-order mixing interaction in- 
duced by ERGE. It is highly expected that the functional form 
and the strength of the entanglement vertex are determined 
in future by these theoretical approaches. In Ref. | |66l] . we 
assumed the following ^ dependence of G by respecting the 
chiral symmetry, P symmetry, C symmetry and the extended 



Figure|9] shows phase diagrams of the EPNJL model. Panel 
(a) presents the phase diagram in d-T plane at 69 = 0, and 
panel (b) shows the phase diagram in 0-T plane at 6* = 0. In 
both diagrams we see that curves of the chiral and deconfine- 
ment transitions almost coincide. 

FigurefTOlshows T dependence of |?7/77o| and \t\}\Q — -k and 
6? = 7r/3. The entanglement vertex also makes P-restoration 
and C violation almost coincide. However, since Tp is slightly 
larger than Tq, there is a naiTow region of T where both of P 
and C symmetries are violated. 

Here, we set Tq to 270 MeV, the original value extracted 
from the pure-gauge LQCD data. Figure niTa) shows T de- 
pendence of a and ^ calculated with the PNJL model at 6* = 
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Fig. 9: (color online). Phase diagrams in (a) T-0 plane at = and 
(b) T-0 plane at 6* = in the EPNJL model. In panel (a)((b)), the 
vertical solid line represents the first-order phase transition caused by 
P violation (C violation). The dashed and dotted lines stand for the 
chiral and the deconfinement crossover, respectively, although they 
almost coincide. 



to 270 MeV with about 50 % reduction of the relative differ- 
ences, - Td and Tp - Tq. 
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Fig. 11: (color online). Results of the PNJL model with Tb ~ 
270 MeV for T dependence of (a) ct/cto and # at 61 = © = and 
(b) |r//r;o| and at 9 — tt and — tt/S. The solid (dashed) line 
represents a (<?>) in panel (a) and rj {xj}) in panel (b). 
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Fig. 10: (color online) (a) T dependence of \r\lr\a \ and aX Q = -k 
and — 7r/3 in the EPNJL model. The solid (dashed) line repre- 
sents rj (?/)). 



and 6* = 0. The difference of the chiral and deconfinement 
pseudocritical temperatures is — ~ 232 — 215 MeV= 
17 MeV. This difference is small compared with the corre- 
sponding result in the case of To = 212 MeV. As seen in Fig. 
[TTl b). however, for 6 — tt and ~ 7r/3 the transition tem- 
perature of P restoration, Tp = 272 MeV, is larger than that 
of C violation, Tq — 238 MeV. The difference is smaller than 
the corresponding difference in the case of Tq = 212 MeV. 
For the PNJL model, thus, the transition temperatures, T^, 
Td, Tp and Tc, are shifted up by changing Tq from 212 MeV 



Figure[T2|a) shows T dependence of a and <P calculated by 
the EPNJL model with ai = as = 0.2 and Tq 270 MeV 
for the case of 6* = and 69 = 0. The chiral and decon- 
finement transitions almost coincide with each other; namely, 
= 222 MeV and Td = 218 MeV. As seen in Fig. [HJb), 
also for 9 ~ TT and = 7r/3, the transition temperature of 
P restoration, Tp = 249 MeV, almost coincides that of C 
violation, Tc = 242 MeV. For the EPNJL model, thus, the 
transition temperatures, T^, Td, Tp and Tc, are shifted up by 
changing To from 212 MeV to 270 MeV with keeping the rel- 
ative differences, T^ — Td and Tp — Tc, small. 

It is widely believed that, in the realistic world, the chi- 
ral and the deconfinement transitions are cross over at 6* = 
6* = 1 84]. However, in a theoretical point of view, it may 
be interesting to consider both chiral and deconfinement tran- 
sitions are of first order. In fact, in the Hol ographic QCD, 
both of the two transitions are of first order 11851 [8611 . (Re- 
cently, the RW transition is also confirmed in the Holographic 
QCD 1 87]. ) Figure [T3]' a) show the T dependence of chiral 
condensate a and Polyakov-loop ^ at 6* = 6* = when we 
put ai = a-z = 0.28 in the EPNJL model. In this case, the 
chiral restoration and the deconfinement transition are of first 
order and happen at the same time. This result is consistent 
with Fig. 3 of Ref. f68tl. Figure fTSl b) show the T dependence 
of rj condensate and the phase tp of the modified Polyakov- 
loop ^ si 6 = 9 and 6* = 7r/3 with the same parameter set. 
Both the P restoration and the C violation are of first order 
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Fig. 12: (color online). Results of the EPNJL model with To = 
270MeV and qi = Q2 = 0.2 for T dependence of (a) a/ao and <^ 
at 6* = e = and (b) 177/770! and It/-] at e = tt and = 7r/3. The 
solid (dashed) line represents a {$) in panel (a) and 77 {ip) in panel 
(b). 



Fig. 13: (color online). Results of the EPNJL model with a\ = 
a2 — 0.28 for T dependence of (a) ct/cto and <P sA9 — & = Q and 
(b) I77/770I and |'i/>| sXO = n and — tt/S. The solid (dashed) line 
represents a (<f ) in panel (a) and 77 (tp) in panel (b). 



and perfectly coincide. There is only one critical temperature 
Tpc ==187MeV below (above) which the P(C) is violated. In 
this case, the P and C violations exclude each other perfectly. 



IV. SUMMARY 

Using the PNJL model, we have investigated P-restoration 
and C violation at finite 9 and/or finite 0, where 6 is the pa- 
rameter of the so-called (?-vacuum, /i = iOT is the imaginary 
quark-number chemical potential and T is the temperature. 

The P symmetry is spontaneously broken below the criti- 
cal temperature Tp at 6 ~ tt and = 0, while the C sym- 
metry is spontaneously broken above the critical temperature 
Tc at = and = tt. The second-order endpoint of P- 
restoration 6 = ttis connected with the chiral crossover, while 
the first-order endpoint of C violation at = tt is done with 
the crossover deconfinement transition. As a consequence of 
these connections, the chiral and the deconfinement crossover 
in the real world with no 6 and turn out be remnants of the 
endpoint of P-restoration and the triple point of C violation, 
respectively. 

When 6 = TT and = tt (vr/S, Stt/S), at zero T, P symme- 
try is spontaneously broken while C symmetry is conserved. 
As T increases, P symmetry is restored just after C symmetry 
is spontaneously broken, so that either P or C symmetry or 
both the symmetries are spontaneously broken for any T. 

Two-flavor LQCD data show that the chiral and the decon- 



finement crossover almost or exactly coincide at 6 = = 0. 
The coincidence suggests that P-restoration and C violation at 
9 ^ ~ TT also almost or exactly coincide. This suggestion 
is supported by the EPNJL model that reproduces LQCD data 
at 6' = 6> = 0, at finite imaginary chemical potential and at fi- 
nite isospin chemical potential. It is natural to think that such 
singular behaviors at 9 = ~ tt take place simultaneously. 
If so, the coincidence between P restoration and C violation 
may be an origin of the coincidence between the chiral and 
the deconfinement crossover in the real world with no 9 and 
0. Further study along this line is quite interesting. 
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